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Abstract: We study one- and two-soliton solutions of noncommutative Chern-Simons 
theory coupled to a nonrelativistic or a relativistic scalar field. In the nonrelativistic case, 
we find a tower of new stationary time-dependent solutions, all with the same charge 
density, but with increasing energies. The dynamics of these solitons cannot be studied 
using traditional moduli space techniques, but we do find a nontrivial symplectic form on 
the phase space indicating that the moduli space is not flat. In the relativistic case we find 
the metric on the two soliton moduli space. 
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1. Introduction 

Soliton solutions of noncommutative field theory (for an introduction to the subject and 
references see for instance 0, |, |, g) is a fascinating subject. In some cases, theories that 
do not have soliton solutions classically turn out to have them when made noncommutative 

The fact that these solitons have an interpretation as lower dimensional D-branes in 
the context of string-theory || |lj and as quasi-particles in the context of noncommutative 
descriptions of solid state systems [§| makes this phenomenon even more interesting. We 
therefore believe that it is worthwhile to study solitonic configurations and their detailed 
dynamics for a range of noncommutative theories j|, [l0| . 

Noncommutative Chern-Simons theory is particularly interesting, as Susskind has ar- 
gued that it describes the fractional quantum Hall effect j/J, (see also [||, 11, 12, |l3[ [bfl). 
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In Susskind's description the solitons are the quasiparticles that exhibit fractional statis- 
tics and are responsible for the unusual quantum behavior of the QHE system. According 
to Susskind, it is the noncommutativity of the theory that encodes the graininess of the 
two-dimensional charged fluid in a magnetic field. 

In this work we study soliton solutions of noncommutative Chern-Simons theory cou- 
pled to relativistic and nonrelativistic scalar fields 15, [l6[j; these theories are noncommuta- 



tive versions of the Jackiw-Pi model |L7], MS] . We find the shape of two-soliton configuration 



using the exact Seiberg- Witten map 20, 21, 22 1. We also find new stationary time 



dependent single soliton solutions, all with the same ^-function charge density profile, but 
with increasing energies. Finally, in the relativistic case, we are able to find the dynamics 
of the two soliton solution using traditional moduli space techniques. 

We think that our results are novel and striking; however, in contrast to the case of 
pure noncommutative Chern-Simons theory, the details of the physical interpretation, in 
particular the relation to the fractional quantum Hall effect, remain to be sorted out. 

This paper is a thoroughly revised version. We thank P. A. Horvathy, L. Martina and 
P. C. Stichel for illuminating comments that led us to reconsider some of our previously 
reported results; see note added at the end of the paper. 

2. Pure Chern-Simons 

We take our noncommutative space to be defined by the commutation relation [x, y] = —id. 
Noncommutative Chern-Simons theory is defined by the Lagrangian 

L = -nK9e a ^Tr(A a dpA J -jA a ApA^ , (2.1) 

where k is the Chern-Simons level, representing the (inverse) filling fraction in the appli- 
cation to the fractional quantum Hall effect 1 . Introducing the notation X 1 = x % — OeijAj 



we write the Lagrangian as [|23|, 24, 25| 

TTK 



L = -— Tr (-EijX* [X* - i [A Q , X>] J + 204)) • (2-2) 
We may introduce the complex notation 

(x 1 — ix 2 ) , 



20 

K=^= (X l -iX 2 ) , (2.3) 



'20 

in which the Lagrangian can be written |l5| 

L = wrKTr (K ] D t K - KD t K^ - 2ttkTiAq , (2.4) 
where we have also introduced DtK = dtK — i [Aq, K\. The equations of motion are 



K, K ] 



1 - 



D t K = . (2.5) 



1 In what follows we shall need to take k positive to find our solutions. 
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Susskind argues @ that in the presence of a quasiparticle the first equation of (p.5[ ) gets 
modified to 



1 + #><0|, 



(2.6) 



where \q\ < 1. In [26] this system was studied from the point of view of String Theory. The 

and 



shape of the solutions was found using the exact Seiberg-Witten map IS, |2C], gl], 
it was shown that the solution represents a deficit or excess of charge centered around the 
origin. However, this solution does not have any moduli, which we need to study moving 
solutions. The situation improves when we include a scalar field, which has been argued 



to describe density fluctuations of the quantum Hall liquid [15|. In this case the solutions 
have moduli describing the positions of the independent solitons, so we should be able to 
make them move around. 



3. The nonrelativistic scalar field 

The Lagrangian for noncommutative Chern-Simons theory coupled to a nonrelativistic 
scalar field can be written as [15, 17] 



L = L CS + 27r0Tr (iDtMt - -^A^(A0) f + H^) 2 ) , (3.1) 
\ 2m J 

where m is a mass parameter of dimension (length) -1 and A is a parameter of dimension 
length. This theory has a BPS bound in the case where 2m\n = 1. Assuming this relation, 
rescaling the time variable t — * mt and introducing the complex notation 

D = - iD 2 ) , (3.2) 

we may write the action in operator notation as 2 

L = otkTt ^K ] D t K - KD t K^ } - 2vrKTrylo 

+27r0Tr |iA# f - (^(^ + D<t>{D<jrf) + ^(^) 2 } , (3-3) 

where D, D and A are defined as 

D t (f> = <p - iA Q cf), 

D<j> = (Kef) - 4>c) , (3.4) 
D<j> = - (K^4> - (pc 1 ^ . 

The action ([T^) is gauge invariant under the transformations 

4>^U4> , K — > UKtf , A -> U(idt + Ao)tf , (3.5) 



2 In fact, a further rescaling is possible which leaves an action with the only free parameter an overall 
factor of re; t — > 8t, (f> — + \J k/6</>, Aq — > Ao/(rnO). 
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where U is an arbitrary unitary operator. The action is also invariant under translations 
by a constant z: 



%D t 



4> -» cpU(z) , K ^ K + z , A -» A ; U{z) = e zc 2c . 
Extremizing the action ( |3.3| ) gives the equations of motion 

, 

l-[K,rf] = -#t , 
{k,^ 1 "} -2(f>c(f) j = 2KiD t K . 
In the time-independent case one can show that any solution to the BPS equations 



(3.6) 



(3.7) 



A 



K, K j 

A = - 

is also a solution to the full equations of motion. 



(3.8) 



2k 



3.1 The one soliton case 
3.1.1 The solution 

The static one-soliton solution discussed in HTBtl looks like 



|o)<4 

K = s|0)(0| +SxcS\, 
1 

29 [ 



(3. 



A 



o 



where we have introduced the shift operators S n = YliLo N + n )(i\ an d the usual coherent 
state \z) = U(z)\0). This solution can be found from the simple soliton at z = by a 



combination of the translation (3J3) and a gauge transformation with a unitary operator 



T(z) 



zK'f-zK 



(3.10) 



3.1.2 The shape of the one soliton solution 

To find the shape of the soliton solution we have to extract it from the operators X\ and 
X2, or equivalently from K = —k= (X 1 — iX 2 ). This we do by using the exact Seiberg- 
Witten map [19, ^C], 21, 22 1, which implies that the Fourier transform of the density is 
given by 



p(k) = Tr (V ifcaXa ) = Tr (e* 



(3.11) 
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where we have defined k = y |(/si — ik?*). We then have 

where x a are ordinary commuting coordinates. This form of the density is very natural 
if one thinks of the coordinates of the individual electrons as the eigenvalues of the X 
matrices, since then p = Yli=i e _lk ' Xi is exactly the Fourier transform of the density p = 

To evaluate the Fourier moments we first notice that K is composed of two pieces 
which commute with each other. This makes it possible to write 



Tr e 



where we have introduced the notation tr for traces over a finite set of states (in this case 
over the set |0)) and Tr' for an infinite trace with some states omitted (in this case |0) is 
omitted). These traces can be evaluated separately. The first one is 

(0| f e~ it ^ g+k ^ ^ ^) |0) = e ~ i( ^ z+k2 ) (3-14) 



and the second one is 

oo oo 
n=l n=l 

Noticing that S\ will shift the state \n) to \n — 1) and that 1 — S\S\ is the projection 
operator on the space spanned by |0), we get 

This trace we perform using coherent states 

Tr (V^+fcct)^ = f f^ {wle -*Ckc+kcr )H 



IT 



^( W \e-^e- ihc \w) e"¥ (3.17) 

7T 

^ 2w c -»(fcw+fcw) c -# 
7T 

= y^ (2) (k) . 

Thus we find the result 

Tr (e-i&K+ktf)} = ^-5^ (k) + e -*(**+**> . (3.18) 

Fourier transforming back we get 

p(x) = — + (5(xi - V2eile(z))5(x2 - V281m(z)) , (3.19) 
2ttu 

which is clearly to be interpreted as a constant background density with a delta function 
type deformation 3 located at \^26z. 

3 Note that the sign of the deformation is opposite to that of the quasiparticles studied by Susskind |Q . 
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3.1.3 More stationary solutions 

The configuration ( |3.9[ ) turns out to be a member of a huge family of stationary solitons 
with a time independent density of the form ( |3.19 ). To see this consider an Ansatz 



4> = ,/tjIoXi* 



Ao = -^|0>(0 



K = S x cS\. 



Inserting (3.20) into the equations of motion (|3.7|) we get 



<v#> 

(V>l<#) 
idM) 



1, 
c'c 



where t = g. 

The general solution of ( |3.23 ) reads 



^ b n e inT . 



n=0 



where 6 n are time-independent and satisfy 



E 

n=0 



)V. 



b n b n +i _ g 



n=0 



If we introduce the generating function 



n=0 



then ( 3.25 ) and ( 3.26| ) can be rewritten as 



7T 



- / d 2 ( G(C)G(C)e-« = 1, 



d 2 CCG(C)G(C)e-« = 0. 



(3.20) 



(3.21) 
(3.22) 
(3.23) 



(3.24) 



(3.25) 
(3.26) 



(3.27) 
(3.28) 



Any function G which satisfies these relations yields a solution to (| 

We may think of the states \ip n ) = e mT \n) as a basis and the combination (3.24) as the 
expression for a general state in this basis. The energy of a configuration can be computed 
from 



E = vrTr {d<P {D(j)) ] + D<p (D<f>) ] 



(3.29) 
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which for our particular ansatz becomes 



E 



TTK 



CC' 



+ 



3 t<#) - 



2-7TK 



The "basis" states thus have energy 



2ttk 



-n, 



(3.30) 



(3.31) 



and a general state has energy 



E 



2lTK 



E| On I 
(n- 1)! 

n=l 



2k. 



A |CG(C)| 2 e-« 



(3.32) 



which can be arbitrary. 

More general solutions can be found by translating these solitons to arbitrary z using 
the transformations (|3.6| ); a simpler form of the translated solution can be found by com- 
bining the translation with the gauge transformation ( 3. 10) ). The naive energy density is 
not gauge invariant; the physical energy density is found by performing a Seiberg-Witten 
map ( |3.11| ). We find that the physical energy density is translated and is concentrated at 
the position of the soliton. 

The basis states \ip n ) can be distinguished from more general configurations in a gauge 
invariant way. If we consider invariants of the form 

/[/]=Tr(0V/( C)C t)) , (3.33) 
then for generic functions /, we find that /[/] is time independent only for the basis states. 



The original soliton (3J3) is the only solution to the BPS equations (|3.8| ) in the family 
constructed here. 



3.1.4 Moving one-soliton solutions 

To make the solitons move one may use the "exotic" Galilean invariance of the theory. 
The infinitesimal version of this tranformation was found in [^7]] . Defining an infinitesimal 
complex paramteter v = -^={v\ — W2) and writing everything in the operator notation, the 
transformations look like 



5K 
SA 



vc 



vc, 

vc' — vc, K 



4 



+ iU(p ( vc^ + re 
+ vt , 



(3.34) 



vc 



vc, A 



+ i 



This transformation can be integrated to finite values of the parameter v. The transfor- 
mation becomes 



<p -> Uct>V ] , 
K -> U{K + vt)U ] , 

A ^U (a + i( K v(K ] - c ] ) - v{K - c))) U ] , 



(3.35) 
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which can be checked to leave the equations of motion invariant. Here U is the unitary 
translation operator which translates a state by vt 



U(vt) = e^-^ 1 , (3.36) 

and V is also a unitary translation operator, but it translates by v(t — iO). It is given by 

V(vt) = e v(t-iO)J-v(t+ie)c _ ( 3 _ 37 ) 



The transformations ( 3.35| ) close under group multiplication up to a U(l) gauge transfor- 
mation: 

?7(uit)Z7(^t)0y+(«2t)V^(«it) = e¥ 2 ^-^)u{{ Vl + v 2 )t)(j)V\{v 1 + v 2 )t) ; (3.38) 

as the U(l) transformation is constant and commutative, it does not act on K and Aq. 

The form of the Galilean tranformation can be improved by acting with an additional 
gauge transformation ( |3.5[ ). In particular, choosing the unitary operator W(vt), one can 
write the action of the Galilean transformation as 

<j> - 4»v ] , 

K -► K + vt , (3.39) 
A ^ A + i (vK ] - vK 



Note the similarity to the translations ( |3.6[ ). 

When K is time-independent, the shift in Aq can then be eliminated by a further gauge 
transformation with 

T = e tvK *- WK , (3.40) 
giving us the final galilean transformation 

<f> -> T(f>V ] , 

K -> TKT ] + ut , (3.41) 
A -> T^oTt ; 

we emphasize that the simple homogeneous transformation of ^4o requires time-independent 
K, 

As an example, the zq = static soliton 

K |o)(o| , 



if = S lC S\ , (3.42) 



= -^|0)(0| 
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can be boosted using (CT) , as K is time-independent. Here, the operator T is block 
diagonal and acts trivially on |0). Also, one may show that TKT^ = K-vt [K, Kl] . Thus 
we find the boosted operators 

* = yf|o>Mf-w)| , 

K = vt\0}(0\+S lC S\ , (3.43) 

A) = -^|o)<o| , 

where \v(t — iff)) is the coherent state at z = v(t — iff). Other solutions can be generated 
by boosting more general stationary solutions, e.g., solutions centered at some point zq or 
the new stationary solutions described above. 



3.2 The two soliton case 
3.2.1 The solution 

After performing the gauge transformation ( |3.5| ) with U = | )( 1 1 + 1 1 )(0 | + the two 



soliton solution found in [15] can be written as 



(A|0)(+| + 5H><-|) (3.44) 



K = z(C\0)(l\ + -\l)(0\)+S 2 cSl, 



where |±) are the normalized states 

1 



= (l*>±l-z», \z) = *r- c |o). 

/ 2 (l± e - 2 N 2 ) 

The BPS equations now tell us that \A\ 2 + \B\ 2 = 1, which is solved by A = sin (a) and 



B = cos (a), that C = cot (a) y coth \ z\ , and that a has to satisfy the equation 

\z\ 2 {c 2 -^j=B 2 -A 2 . (3.45) 

This is a cubic equation in tana, so it can be solved explicitly. The solution is not par- 
ticularly illuminating, so we explicitly give only some limiting properties and use the full 
result in the calculations. We find that 



e -2\z\ 2 

c = 1- 

4 M - 

C = \z\ (l + \ N 4 - H \ z f + O ( \z\ 12 ) ) 1*1 « 1 . (3.46) 



- + 0(e-^) |*|»1 
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3.2.2 The shape of the two soliton solution 

Using the same method as in the one soliton case we can calculate the shape of the two 
soliton solution. The Fourier moments of the density are given by 

p(k) = Tr ( e -i(kK+kKi)\ ( 3 _ 47 ) 



which again splits into two pieces 

rft ( e -i(kK+kK^)\ = tr / e -i(D|0>(l|+5|l>(0|)\ + Tr / (-i{kS 2 cSi+kS 2 Jsi)\ ^M) 



where we have introduced D = kzC + |f and where tr now means the trace over the 
subspace spanned by |0), |1) and Tr' means the trace over the orthogonal complement of 
this space. Using the methods developed for the one soliton case we can evaluate this as 

Tr' ( e -i(kS2cSl+kS 2 Jsl)\ _ Tr ^ e -i(kc+kci)^ ^ ( 3 _ 49 ) 

which again leads to a constant background density, and 

tr ( e -^|0>(H+5|i>(0|)\ = 2cos{m f (3 . 50) 



which gives us the perturbation part of the density 

P P (x) = J ^2cos(|L>|)e^ . (3.51) 

To do this integral we choose integration variables k that simplify \D\: we change coordi- 
nates from k to q 

' gi +^, (3-52) 



2\z\ \X+ A 



A ± = in ( 1 ± C ) . (3.53) 



where we have defined 

\ , = 

2 \C 

In these coordinates D is simply 

\D\ 2 = ql+ql. (3.54) 
Using the asymptotic expressions for C ( 3.46| ) we see that 



| \z\ « 1 
\z\ \z\ > 1 



1 



\z\ « 1 



A = < 2 (3.55) 
I \z\ > 1 V ' 



The Jacobian is 2 8\ + \ so we nave 



^> = e6rJ^ co i^) e "" v " <3 ' 56) 
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where we have introduced a new variable y defined by k • x = q • y, which implies 



(y 1,2/2) 



1 



4 Id A. 



-[x\{z + z) + ix 2 (z - z) 



1 



4 U| A_ 



[x 2 (z + z) — ixi(z — z)] ) . (3.57) 



Note that a circle of radius 1 in the y coordinate system when transformed to x coordinates 
looks like an ellipse with the major axis oriented in the direction of z and of length ^/29\ + 



and a minor axis of length V2#A_. From ( 3.46 ) we also see that for small \z\ the ellipse 



becomes a circle of radius \j §, whereas for large \z\, A_ goes exponentially fast to zero so 

±V29\z\. 



the ellipse is rapidly squeezed into a line between ±V2#Ah 
We perform the angular integral and find 



Pp(y) 



where y is the absolute value of the vector y. This integral is divergent and needs to be 
regulated. This we do by adding a small real part to the argument in the exponential. 
Thus, with a = i + e we have 



Pp{y) = liir A 



1 



ex , a. 



-Tie 



dq 
2^ 



qe aq J (qy) = lim 



1 



Ke- 



2vr6lA + A_ ( y2 + a 



2^1 



(3.59) 



which is completely well defined and smooth for all non-zero e. Letting e — > the function 

is zero for y > 1 and for y < 1 it looks like p p = — 27T e\ \ ~ — r • I n figure || we see a 

cross section of the full density drawn in the y-plane where the border of the ellipse is at 
\y\ = 1 and for a small non-zero value of e. The background density is normalized to one. 



—\ 1 

1.5 2 
lyl 



Figure 1: The full density 



Notice that the full density inside the ellipse is always less than the background density. 
This means that we have taken particles from the inside of the ellipse and piled them up 



at the border of the ellipse. Notice also that there is always a small region on the inside 
of the border of the ellipse where the full density will be negative. This seems difficult to 
reconcile with the interpretation that the theory describes an electron gas since then there 
could not be less than zero electrons at any point. 

As a final exercise we show how two one-soliton solution can be found in the \z\ — > oo 
limit of the two-soliton solution. We want to show that 

lim p p (x) = <5 (2) (x - VWz) + 5 {2) (x + VWz) (3.60) 

\z\— >oo V / V / 

To do this we have to use the regularized version of p p using the non-zero e which we will 
take to zero only at the end of the calculation. Thus we would like to show that for any 
real function f{xi,X2) 



f(V26z) + f(-V26z) = lim lim f d 2 xf(x) Pp (x) 



= lim n^ im 9-71-/9 X X Ue / d V(x) ^r(3-61) 

Changing integration variables from x to y, not forgetting the Jacobian, we get 

lim lim -Tie [ d 2 y f(V20X +yi , V26X-y 2 ) — / (3.62) 

But in the limit when \z\ goes to infinity, A_ goes to zero exponentially fast. This means 
that the y 2 dependence of the function / is scaled out and the integral becomes 



1 f f+oo ^ 
lim lim -Tie / dy ± /(v^A+yi, 0) / dy 2 T 

e^0| 2 |-Kx>7T J 7_ 00 ( y 2 + (j + e )2)2 

= lim lim -Tie [ d Vl f(VWX +yi ,0) * (3.63) 

Taking the real part and performing the e — ► limit we recognize a representation of the 
delta function 

limfte 2 2t = lim ^ = 2vr ( 5(y 1 2 - 1) = tt%i - 1) + nS( yi + 1)(3.64) 

and substituting this back we get 

/(v^A + ,0) + /(-v^A+,0) (3.65) 
and since for large \z\ we have A + « \z\, this is what we wanted to show. 
3.2.3 The dynamics 

To study the dynamics on the two soliton moduli space one can try to make the moduli 
slowly time dependent and find the effective action in the moduli. Doing this we find 



2 

2 



i-h jdl^C- ^) |:.|- - z~) . (3.G0) 
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or, using z = re 1 ^, 



-2~k / <// I C z - ] ,-■ 



(3.67) 



Notice that there is no time derivative of r and that the function in front of 4> depends only 
on r. It can be written as (4A+A_) 2 which is identical to the "angular momentum" of the 



two soliton configuration defined in |15| . Using the asymptotic expansions for C ( 3,46 ) we 
find that it goes as (1 — 4r 4 ) for small r and for large r it goes exponentially to zero. 

The effective action is first order in time derivatives and thus comes in "Hamiltonian" 
form (with the Hamiltonian equal to zero). It defines a connection 1-form on the moduli 

c 7 ) 



space A 



C 2 



(zdz — zdz) whose curvature 2-form is the symplectic form 



on the moduli space seen as a phase space. This may be taken as the starting point for 
quantization of the system but since the Hamiltonian is zero, there is no non-trivial time 
evolution. We have plotted the symplectic form in figure At small \z\ it goes to zero as 




Figure 2: The angular momentum and the symplectic 2-form 



0.6) there is a 



16|z| and at large \z\ it goes to zero exponentially . In between (at 
maximum. 

One could speculate that, as in [28], a potential for the moduli could be generated. 
This would then define a nontrivial Hamiltonian which would give non-trivial dynamics on 
the moduli space. 

As found in |2£ 



29, 3C] in the commutative Chern-Simons case, the symplectic form 



defined above equals the Kahler form of the metric on the moduli space of vortex solutions 
to the Maxwell-Higgs model. If this is true also in the non-commutative model, the sym- 
plectic form we just computed should have the properties of a Kahler metric. This seems 
to be true for small \z\ where the apparent conical singularity is removed after dividing by 
the Z2 action that arises because the voritces are identical particles. At large \z\ however, 



'Since the form degenerates at and 00, the phase space is singular at these points. 
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the metric is zero, which seems to go against the intuition that the metric should be flat 
for well separated vortices (they should not interact at large distances). It also seems to 



disagree with explicit calculations of the metric in the Maxwell-Higgs case [31]. However, 
our result is similar to what we find in the relativistic case in the next section, where we 
again find a nontrivial metric between vortices even in the large \z\ limit. As in |3(]] one 
should be able to use our expression for the symplectic two form to calculate the phase 
space volume or rather the reduction of the phase space volume due to the presence of 
other vortices. From this one could find the statistics of the vortices. 

There have been other attempts to find the dynamics of first order system of vortices, 
e.g., [ 33 1 . However, without additional assumptions the methods of do not yield unique 
solutions, and we were not able to put them to fruitful use in our case. 

4. The relativistic scalar field 

The Lagrangian of a noncommutative Chern-Simons field interacting with a relativistic 
scalar field can be written |l5| (using the objects defined in (|3.4|)) as 



L rcl = ivrKTr ^K j D t K - KD t K j } - 2ttkTtA + 

+2vr0Tr {a^A^ - \ (D4>{D^ + D4>{D4>)^ - V(^) } (1.1) 
leading to the equations of motion 



1 



DM--{D,D}ct> + V l 



, 



KiDtK 



2(f>cq 



(4.2) 



(4.3) 



1 - [K,K ] ] = i- 

K 

If the potential has the form 



A 



■(A 



v(o = i,ac-u 2 f , 



(4.4) 



(4.5) 



where u is a (real) constant, then the model has a Bogomolny bound and BPS equations: 

D<j) = , (4.6) 



At _ , ( 2 



(4.7) 



1 - [K, K^] 



i— 

K L 



A<; 



■(A 



(4.8) 



There is also a second set of BPS equations in this model, with D(j) = and A 



u 



! ) (f), but these equations do not seem to have localized solutions as there are 



no normalizable eigenvectors of the creation operator c 
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It is instructive to see how ( fO^JI^fO]) follow from ( |4^j47| , |48|) . Substituting (| 
into (fOI) we get 



1 - [K, K^) 



29 



-u 2 ) 



From D(j) = and the definitions of D, D it follows that 

-\{D,D}<p = ^[D,D]4> = ](l- <t> = ~^ (V " u 2 ) #fy 

Computing the D t derivative of (|4.7|) , we get 

D t D t <f> = -- (u 2 D t <P - D t <f><t>s<t> - <p(Dt<P)U - 4>4> ] D t( ^ ; 

hence 

AA^ = — (# L «jV 



and thus 

u Ki 

For V of the form (fQ|), this is Q. 
To derive (|4.3|) , we write 



-V' 



= AW) = -(D t K)^ - D t (f> + D t ^ 



which gives 



*D t K = c(A^) t - {D t ^K = — | r, 



multiplying with on the left, we get 



i> f ikD+K = 6cd) ] 



Noting that D(f> = =>- 



this gives 



fct 



{ir,# t }-2^ t 



(4.9) 



(4.10) 



(4.11) 



Here we see that, unlike in the non-relativistic case, a solution to the BPS equations imply 
a solution to the full equations of motion and not just to the time independent equations 
of motion. This is almost true. There could also exist BPS solutions where the operator 



i K D t K-{K, # T } + 2<Pccp\ 

is non-zero, but lies in the subspace orthogonal to the operator 4>($ . These "solutions" 
would then be solutions to the BPS equations only. 
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4.1 The one soliton case 



4.1.1 The static solution 



The static, one soliton solution is given by the Ansatz [15] 



K = z|0>(0| +S lC Sl, <j> = \\0)(z\, 



An 



I o ><0 | 



where the parameters A and b are real and are determined in terms of the parameters k, 9, u 



in the Lagrangian. It follows from the equations of motion (4.2.4.3,4.4) that b = and 
A 2 1 



1± A 1 




A 2 1 / , I 6k 2 \ 
OT ^ = 6 1:t V 1+ ft? • 



The first two values correspond to configurations that are also solutions in the Bogomol'nyi 
limit. Note that they exist only for 9 > 2k 2 /u 4 , whereas the remaining two solutions are 
well defined for any nonvanishing noncommutativity parameter. 

The K part of the solution has the same form as in the nonrelativistic case and the 
shapes of relativistic and nonrelativistic solitons coincide. 

4.1.2 More stationary solutions 

Encouraged by the results from the nonrelativistic model, one may try to find stationary 
solutions with the same charge density as the BPS solutions. In this case however, one 



would have to use noncommutative Lorenz boosts [35] instead of a Galilean boosts to make 
the solitons move. We have not worked out the details of this but we feel confident that 
the solutions can be found. 

4.2 The two soliton case 
4.2.1 The solution 

As in the nonrelativistic case the two soliton solution can be written in the form [ff5ll 



<f> = A (sin a | ) ( + | + cos a | 1 ) ( — | ) , 
K = z(c\0)(l\ + ±\l)(0\)+S 2 cSl 



(4.12) 



and the BPS equations (4^6,4/7,4^8) give 



1 - |z| 2 (C 2 - C~ 2 

i + M 2 rc 2 -c- 2 ) 



19\ . n / 2 \2 ■ 2 \ 

■ — — sm a [u — A sm a] 

20A 2 n / 2 \ 2 2 \ 

— 2— cos a [u — A cos a) 
cot 2 a = C 2 tanh |z| 2 



(4.13) 



Unlike the nonrelativistic case, it is not possible to solve ( 4,13j ) in a closed form. It is how- 
ever straightforward to find approximate solutions for small and large soliton separations. 
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We find 



C 



cot a 



\z\ 2 + 



^-\)\*\<+o(\*n 



\2 
A 



2Xt 



u 2 - 2A 2 



\z\ A + 0{\z\*), 



4^2 



for Id <C 1 and 



C 

cot a 
A 2 
A 2 



flA 2 ^ 2 - Xl) 



4«2| z |2_ eA 2 o(2u 2_ A 2 o) 

4k 2 |z| 2 | 



4 K 2 |z| 2 -0A 2 c > 2 -A 2 c 

A 2 + O(e- 4 N 2 ), 



e -2|,| 2 + (V 4 ^ 



it 



1± A/1 



r 



(4.14) 
(4.15) 
(4.16) 

(4.17) 

(4.18) 

(4.19) 
(4.20) 

(4.21) 



for \z\ S> 1. 

Note that — modulo different dependence of the coefficients on z — the form of K 
in relativistic and nonrelativistic cases coincide and the leading terms in the C expansions 
for small and large \z\ are the same. Therefore the analysis of the shape of the two soliton 
solution performed in section 3.2.2 can be repeated for the relativistic solution virtually 
without any changes. 

4.2.2 The dynamics 

To understand the dynamics of slowly moving solitons we use Manton's method of moduli 
space approximation (3^, £57| ( see a l so [H> HH)- The effective action, obtained from (4.1) by 
inserting the ansatz ( 4.12; ) with time dependent moduli z has the form 



(4.22) 



with 



(5 2 A) 2 + A 2 



(d z a) 



(tanh 2 \z\ 2 sin 2 a + coth 2 \ z\ 2 cos 2 a) 



9zz = 2|cU| 2 + 2A 2 



i2 \ ^ \ / 2 i 1 2 2 2ii2 2 A 

\9z a \ t- (tanh |z| sin a + coth |z| cos a) 



+ A 2 Hzl 2 + tanh \ z\ 2 sin 2 a + coth Izl 2 cos 2 a] 



(4.23) 



Using ( 4.13 ) it is straightforward to express d z \ and d z a in terms of A, a and C, but the 
result is complicated and not particularly illuminating and we refrain from giving it here. 



We also not analyze the dynamics that follows from (4.22) in detail and discuss only the 
limiting cases of small and large \z\. 

For \z\ <C 1 the effective action takes the form 



S cS = 2-kOXI I dt r 2 



2u 2 

3 (2r 2 + r 2 ^ 2 ) H r 2 ip 



+ (rV,r 8 ) 



(4.24) 



where z = r and we have dropped a total derivative term. In the coordinates 

R = r 2 , ■&= V2<f , 

this action reads 

S eS = 37r6\ 2 J dt (^R 2 + RH 2 + R 2 ^j 
and leads to the equation of motion 

R- R (u 2 -^-)i> = o 



(4.25) 



d_ 

dt 



R 2 + 



3 k 

s/2u 2 
3k 



(4.26) 



R 2 



V2u 2 
3k 



where c\ is an integration constant. 

The simplest solution of fl4,26| ) correspond to solitons circulating around each other 
with fixed angular velocity, 



R = const , i9 



2 v / 2u 2 
3k 



r = const , ip 



2n2 
3k 



Another simple solution is obtained for c\ = 0, 

,2 



4 2-2 

r = c 2 sin 



3k 



it - 1 ) 



u 

3k ' 



and describes head-on scattering of the solitons. 

In the general case of c\ ^ it is possible (neglecting a term j^f 4 in comparison with 
c 2 r -4 , which is consistent with the approximation made in deriving fl4.24j )) to write down 
the solution in the form 



4(t-t ) 2 + ^ 



.2 ' 



<P 



U Z 1 
— + 



ci4 



3k y^2 c\ + 4 (t - t ) 2 



In the region \z\ 3> 1 one gets 



S cS = 2tt9XI I dt r 2 (r 2 + 2rV) + O (e^^) 



(4.27) 
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In the variables 



R = r 2 , i? = 2V2tp , 



(4.27) takes a form of a free particle action. Consequently the general solution of the 
equations of motion that follow from ( |4 .27 ) can be presented in the form 



Cot + Cl 



with (complex) constants Co an d Ci- 

We thus see that the naive moduli space has conical singularities in both the short and 
long distance limit. This is so even for identical particles where we have to quotient the 
moduli space by a Z2 action, which for the original GMS solitons removed the apparent 
conical singularity at the center of the moduli space |38j. However, the result agrees 
with what we found in the non-relativistic case where the phase space is singular both in 
the long and short distance limit. 



5. Conclusions 

We have studied one and two soliton solutions in noncommutative Chern-Simons theory 
coupled to nonrelativistic and relativistic scalar fields. We studied the shape of the solitons 
using the exact Seiberg-Witten map. We have found new stationary solutions with time- 
independent charge and energy densities. The most general state could be written as a 
certain superposition of basis states. We have also integrated the infinitesimal Galilean 
transformations of 27], and used them to generate moving solitons. 



For the two soliton solutions we tried to study the dynamics using the moduli space 
approximation. The direct approach is possible only in the relativistic case because in the 
nonrelativistic case the action is linear in time derivatives. In the relativistic case, we found 
explicit expressions for the metric in the limits where the solitons are far apart and when 
they are close. In these regimes the metric is flat with a conical singularity at the center. 
However, the conical singularity is different in the two limits. This makes the interactions 
"long range" and the solitons affect each other at large distances. This awkward long range 
interaction could be avoided, but only at the price of an equally awkward redefinition of 
the physical variables. In the nonrelativistic case, we could not find a metric, but found a 
symplectic form that has singularities that resemble the singularities of the metric in the 
relativistic case. 

Noncommutative theories are known to link UV and IR behavior; the long range 
interaction that we find could be related to this, and could be probed by studying the 
theory on a cylinder or torus. The physical meaning of the model should be elucidated by 
coupling to a physical external electromagnetic field. 

The new types of solutions require further study. We should investigate the meaning 
of the gauge-invariant observables that distinguish the basis states with quantized energy 
from the general mixed states. We should also determine if such solutions arise in other 
theories, e.g., the noncommutative Maxwell-Higgs system. 
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Note added: 



An earlier version of this paper found only a special subfamily of the solutions that we 
discovered here, and were led (falsely) to conclude that no Galilean transformations are 
possible, and that moving solitons travel with quantized velocities. We thank P. A. Hor- 



vathy, L. Martina and P. C. Stichel for directing our attention to [27] and enlightening email 
exchange on the subject. We have been informed by the authors that they have also written 
a new paper on related topics [S3]. We also thank D. Bak for useful communications. 
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